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Robust stationary entanglement of two coupled qubits in 
independent environments 
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Abstract. The dissipative dynamics of two interacting qubits coupled to independent reservoirs at nonzero 
temperatures is investigated, paying special attention to the entanglement evolution. The counter-rotating 
terms in the qubit-qubit interaction give rise to stationary entanglement, traceable back to the ground state 
structure. The robustness of this entanglement against thermal noise is thoroughly analyzed, establishing 
that it can be detected at reasonable experimental temperatures. Some effects linked to a possible reservoir 
asymmetry are brought to light. 

PACS. 42.50.Lc Quantum fluctuations, quantum noise, and quantum jumps - 03.65. Yz Decoherence; open 
systems; quantum statistical methods - 03.65.Ud Entanglement and nonlocality - 03.67.Mn Entanglement 
measures, witnesses, and other characterizations 



1 Introduction 

The entanglement dynamics of open quantum systems has 
attracted considerable attention over recent years. Quan- 
tum correlations are indeed central in a variety of con- 
texts, both because of their fundamental properties and 
because they are key ingredients in the fields of quantum 
communication, computation and information processing. 

Superconducting devices are among the best candi- 
dates for quantum state engineering. The nonlinearity of 
Josephson junctions can be used to mimic two-level quan- 
tum systems (qubits), with the advantages that the in- 
herently low dissipation of superconductors make long co- 
herence times possible, while integrated circuit technology 
allows for scaling to large and complex systems [IJ[2j[3]. 

Recent experiments have successfully demonstrated the 
existence of quantum coherent oscillations in individual [H 
\MMM\ and coupled [9j[l0l[IU[l2l[l3] Josephson qubits. 
Yet many aspects of decoherence and dissipation in these 
systems still remain to be understood. A current challenge 
is to discriminate among different sources of decoherence. 
This calls for finding appropriate representations of the 
interaction between these devices and the external envi- 
ronment, so as to single out the conditions (quite often 
model-dependent) under which a coherent behavior is still 
visible. 

Within this scenario, the simple model of two cou- 
pled qubits in thermal contact with different baths is more 
than an academic curiosity, for it can describe quite differ- 
ent experimental setups as, for instance, two inductively 
(capacitivcly) coupled flux (charge) qubits, provided that 



their spatial separation is large. It also allows one to in- 
vestigate the mechanisms that determine the degradation 
of coherent Rabi oscillations in Josephson phase qubits 
stemming from spurious microwave resonators. These res- 
onators, as shown by Martinis and coworkers |14j . arise 
from changes in the junction critical current produced by 
two-level states in the tunnel barrier. One bath describes 
then the electromagnetic environment of the phase qubit, 
while the other one mimics the phononic environment cou- 
pled to the junction microwave resonator. Finally, such 
an approach is of relevance in other frameworks, such as 
cavity QED [T51l 16 [11711 18| . quantum dots, and spin sys- 
tems [TW20] . 

Quite recently, we solved the zero-temperature dynam- 
ics of a model that includes the usually omitted counter- 
rotating terms in the qubit-qubit interaction. Besides con- 
firming many features previously reported in the litera- 
ture [^[1I^[M1CIM^^ our anal- 
ysis showed that these terms produce long-time stationary 
entanglement |35U36[|3"T] . The next natural task is thus to 
look for effects stemming from nonzero reservoir temper- 
atures. 

In this vein, we work out here the solution of the mas- 
ter equation at generic reservoir temperatures, establish- 
ing that the effects of the counter-rotating terms turn out 
to be very robust against thermal noise, since they prevail 
at temperatures at which Josephson qubits work in actual 
experiments. We put forward the role of the qubit-qubit 
coupling constant A in the creation of stationary entan- 
glement at finite temperatures: the optimum case can be 
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reached for values of A of the order of the qubit Bohr 
frequencies, a condition easily met in practice. We finally 
show how these results are modified by possible asymme- 
tries between the reservoirs. 

The paper is structured as follows. The decay model, 
the dynamics at zero temperature and the comparison 
with the results for a single common reservoir are pre- 
sented and discussed in Sec. II. The solution at generic 
reservoir temperatures and the analysis of the entangle- 
ment dynamics is given in Sec. Ill for different initial con- 
ditions. We also investigate the effects of different reser- 
voir temperatures and qubit-reservoir coupling strengths. 
Finally, our conclusions are summarized in Sec. IV. 



2 The decay model and the dynamics at zero 
temperature 



2.1 The model 



dynamics of the qubit-qubit system turns out to be: 
g(t) = -i[H s ,g(t)} 

+c I (\a){b\g(t)\b){a\-\{\b){b\,g(t)}) 
+c n (\a){c\g(t)\c){a\-^{\c){c\,g(t)}) 
+c I (\b)(d\g(t)\d)(b\-±{\d)(d\,g(t)}) 

+ ClI (\c){d\ Q(t)\d){c\-±{\d){d\,Q(t)}) 

+cj(\b)(a\g(t)\a)(b\~±{\a)(a\,g(t)}) 

+c n (\c){a\g(t)\a){c\-^{\a){a\,g(t)}) 

+c I (\d){b\g(t)\b){d\-±{\b){b\,g(t)}) 

+c n (\d)(c\ g(t)\c)(d\-±{\c)(c\,g(t)}) 

+c cr j(\a)(b\ g(t)\d)(c\ + \c)(d\ g(t)\b)(a\) 
+c cr ,n(\a){c\ g(t)\d)(b\ + \b)(d\ g(t)\c)(a\) 
+c cr j(\d)(c\ e (t) |o><6| + |6><a| |c><d|) 
+c clJI (\d)(b\g(t)\a)(c\ + \c)(a\g(t)\b)(d\). 



(2) 



The states appearing in the master equation, namely the 
eigenstates of the Hamiltonian (TTJ) , are given, in the un- 
coupled basis {|00), |11), 1 10) , |01)}, by 



The system under study consists of a pair of coupled 
qubits. Denoting by |0)i (|0)a) the ground state of the 
first (second) qubit and by |l)i ( 1 1) 2 ) the corresponding- 
excited state, the Hamiltonian is given by (with Ti = 1 
throughout) [35]: 



(2) (2) 



(1) 



\a) = cos ^-|00) -sin , 
|6)=cos^|10)-sin^|01) 
| c )=sin^|10)+cos^|01) 
\d) = sin y|00)+ cosy |11), 



(3) 



uji being the Bohr frequency of the ith qubit and A the 
coupling constant, which, in the case of flux qubits with 
flux-flux coupling, is proportional to their mutual induc- 
tance [38] . The Pauli operators are given by o ( + = |l)jj(0|, 

o-W = |0)i<(l| and ai l) = cr® + a®, with i = 1,2. Note 
that ([1} contains counter-rotating interaction terms. 

In reference [35] . a master equation was derived for 
the case in which each qubit interacts with an indepen- 
dent bosonic thermal bath. Indeed, by using the general 
formalism given in reference [39], under the secular ap- 
proximation and the assumption of independent bosonic 
reservoirs, the Markovian master equation governing the 



and correspond respectively to the eigenvalues (ordered 
by increasing energies): 
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where, for L02 > Wi, the parameters 9i and On satisfy the 
relations 

. „ IAI 
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The decay rates c 2 ; and the cross terms c cr ^ (with i = I, II) 
are 



ci = 7/41 cos "77 
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where 



(0) 



(7) 



7i,« = [l + JVi(w 4 )], 

J/(w) being the zero temperature spectral density of the 1- 
th reservoir, and Ni(oji) the number of photons in a mode 
of frequency u>i of the same reservoir. The corresponding 
excitation rates Ci can be obtained by substituting, to n t u 
in the corresponding Ci, the expression: 

'^ /KbT ' = Ji(ui)[Ni(ui)]. 



7i, « = 7i, H e 



(8) 



Since Ni(uJi) — > for 7) — > 0, the latter equation clearly 
shows that the excitation rates vanish at zero tempera- 
ture. We stress that when both reservoir temperatures are 
zero, the excitation rates vanish, which physically trans- 
lates the impossibility of creating excitations in the system 
due to the interaction with the reservoirs. 

2.2 Dynamics at zero temperature and comparison 
with previous literature 

In reference |35| we calculated the two-qubit dynamics 
at zero temperatures. As expected, all the coherences os- 




Fig. 1. Dynamical evolution of the concurrence, at zero tem- 
perature, when the system starts from the states (a) ^/p\01) + 
VT^pllO) and (b) ^/p|00) + ^/T^p\U) as function of time 
and of the dimensionless weight p. 



cillate with an envelope gradually decreasing to zero. As 
for the populations, the state \d) decays exponentially to- 
wards the states |c) and \b), which in turn decay towards 
the ground state \a). From the time evolution of the pop- 
ulations and coherences, the concurrence |40j can be de- 
duced. Here, we focus on two wider classes of initial con- 
ditions. 

In figure l.a the system is initially prepared in the 
one-excitation state ^101) 



v/1 — p\10) and we plot the 
concurrence as a function of time and of the weight p. In 
figure l.b, the system is in a superposition of the states 
|00) and |11), namely y/p\00) + vT^~p|ll). 

For the qubit frequencies we take lo\ — L02 = 5GHz, 
which are typical values for the current superconducting 
technology [TJ03I3]. Since experiments indicate that the 
qubit-qubit coupling constant A can be of the order of 
the single qubit frequencies [I0ll37j . we take A = cji = W2, 
in order to maximize the amount of stationary entangle- 
ment [55] . 

In addition, we consider Ohmic reservoirs with zero- 
temperature spectral densities of the form 



J,(cj) = a 4 w, 



(9) 



and equal qubit-reservoir coupling strengths a\ = ct2- 
This implies that, since "fi,u (711,0) is proportional to the 
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spectral density Ji(wj) [Ji{uJii)\, the ratio between the de- 
cay rates is equal to cu/cj = luu/ujj. Finally, we assume 
a\ = Q2 = 10~ 3 wi , since it reproduces the decay times 
found in recent experiments on superconducting qubits [TJ 

In figure l.a Rabi oscillations due to the coherent ex- 
change of excitations between the two qubits can be ob- 
served at short time. A different yet oscillatory behavior is 
seen in figure l.b, due to the counter- rotating terms in the 
qubit-qubit interaction, as well as the phenomena of en- 
tanglement sudden death [2"Tll2"2"ll2"5] and birth [25] , whose 
origin can be recognized in the different role of coherences 
and populations in appropriately manipulated expressions 
for the concurrence. 

Concerning the long-time behaviour, both figure l.a 
and l.b show that the concurrence reaches a stationary 
value. In reference [26], for example, the presence of this 
entanglement when there is only a common reservoir is 
justified in terms of the very different timcscales in the 
decay of the symmetric state (|01) + \10))/^/2 and the 
antisymmetric one (|01) — |10))/\/2- In fact, this causes 
the strong dependence of the long-time correlations on 
the initial state. 

The situation in our case is quite different. The ori- 
gin of the stationary entanglement can be traced back 
to the structure of the eigenstates of the Hamiltonian 
([1]). The zero-temperature dynamics is such that the sys- 
tem eventually goes to the ground state |a), which is a 
superposition of the zero-excitation state |00) and the 
two-excitation state |11), the latter having a probabil- 
ity amplitude proportional to the coupling constant A, for 
A <C cji,W2 [see equation (J5])]. The net result is a nonzero 
stationary value for the concurrence, irrespective of the 
initial state. 

In principle, this independence from the initial state 
can be exploited to experimentally discriminate between 
the entanglement created by a common reservoir and the 
one due to the structure of the ground state, which is 
a consequence of the presence of counter-rotating terms 
in the system Hamiltonian ([1]). Indeed, if there were no 
counter-rotating terms [that is, a qubit-qubit coupling of 

the type X(a^cr^ + a^a^)], the ground state of the 
system would be simply |00), leading to no stationary en- 
tanglement in the absence of a common reservoir correlat- 
ing the two qubits. 



equations: 

Qaa = —{Cl + Cli)Qaa + Cj Qbb + CllQcc, 
Qbb = ciQaa - (cj + c n )g bb + CnQdd, 
Qcc = CnQ aa - (cj + C U )Qcc + CiQdd, 

Qdd = CnQbb + ciQcc - (c/ + c H )g dd . (10) 

This can be solved by Laplace transforms, which leads us 
to the following algebraic linear equations 

(s + c/ + c n )g aa - a Q bb - c n g cc = Q aa (0), 

-CiQaa + (s + C/ + Cn)Qbb - CnQdd = Qbb(0) , 
-CllQaa + (S + C/ + Ch)Qcc - Cl Qdd = Qcc(0), 

-cnQbb - ciQcc + (s + ci + c H )g dd = Q d d{0), (H) 

where 

r + co 

Qkk = / e~ lst Qkk (t) dt. (12) 
Jo 

After some algebra and transforming back to the time 
domain, we finally obtain the solutions 



3 Dynamics at arbitrary reservoir 
temperatures 

3.1 Equal reservoirs 

Next we consider the nonzero temperature dynamics. The 
solutions for the coherences given in |35j are still valid 
for arbitrary reservoir temperatures. On the other hand, 
from equation ©, it is possible to show that the dynam- 
ics of the populations is governed by the following rate 
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Qaa(t) 



cicn 



(cj + ci)(c n + c n ) 
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Fig. 2. Dynamical evolution of the concurrence, at temper- 
atures (a) Ti = T 2 = 10 mK and (b) 7\ = T 2 = 20 mK, 
when the system starts from the state ^/p\01) + \/l — p\ 10} as 
function of time and of the dimensionless weight p. 



which reproduce, at t = oo, the stationary state at generic 
temperatures and, for c/ = cjj = 0, also the zero-temperatur 
solutions given in Ref. [33] . 



^From equations (|14| it is possible to evaluate the ef- 
fects of nonzero reservoir temperatures on the entangle- 
ment. In the following we assume T\ = T2 = T, and again 
equal coupling strengths between each qubit and the cor- 
responding reservoir. We postpone to the next subsection 
the analysis of the effects of different bath temperatures 
and different coupling strengths. Figures 2 and 3 show 
the time evolution of the concurrence, as a function of 
the time t and the weight p, when the system starts from 
the states ^01) + yT^p|10) and 7p|00) + s/T=p\ll), 
respectively. 

In figure 2 we can see that the oscillations in the con- 
currence, which mirror the coherent periodic exchange of 
energy between the two qubits, are quite robust for short 
times against thermal noise. Indeed, even at the rather 
high temperature of 20 mK, they remain visible for a 
few Rabi periods. The same happens to the short-time 
qubit-qubit coherent oscillations, caused in this case by 
the counter-rotating terms when the system starts from 
the state v /p| 00) + ^Jl — p|ll): again the oscillations are 
clearly visible at reservoir temperatures of 10 and 20 mK 
as shown in figure 3. 

e From both figures 2 and 3 it is apparent that the sta- 
tionary entanglement is less robust against thermal noise 
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Fig. 3. Dynamical evolution of the concurrence, at tempera- 
tures (a) Ti = T 2 = 10 mK and (b) Ti = T 2 = WmK, when 
the system starts from the state ^/p|00) + \f\ — as func- 

tion of time and of the dimensionless weight p. 



the stationary state is very close to a maximally entan- 
gled state when A — > oo, for experimentally meaningful 
temperatures (i.e., of the order of tens of mK) the opti- 
mal value for stationary concurrence corresponds to those 
A up to ten times the qubit Bohr frequencies. Indeed, for 
A — > oo the value of C(oo), which at zero temperature is 
almost unity, decays very rapidly with temperature. 



3.2 DifFerent reservoirs 

Let us conclude our analysis bytaking into consideration 
different reservoir temperatures and qubit-reservoir cou- 
pling strengths. 

Figure 5 shows the stationary concurrence as a func- 
tion of the reservoir temperatures T\ and T2, for A = io\ = 
u>2 and equal qubit-reservoir coupling strengths. Figure 6 
compares, in terms of the corresponding contour plots, this 
situation with the case in which the coupling strengths are 
different. 

For equal coupling strengths, we note an almost sym- 
metric role of the reservoir temperatures in the destruction 
of the stationary entanglement. From figure 6 we note that 
the only effect of different coupling strengths is to erase 
this symmetry, in the sense that the larger the coupling 
strength, the smaller the temperature necessary to repro- 
duce the same damping. This interplay between temper- 
atures and coupling strengths can be easily understood 
if we recall that the two reservoirs give independent con- 
tributions to the decay and excitation rates and the dy- 
namics is influenced only by the numerical values of these 
rates. 




Fig. 4. Stationary concurrence C(oo) as a function of the bath 
temperatures Ti = Tb = T and the coupling constant A. The 
values of A giving robust stationary entanglement are of the 
order of 10 GHz, i.e. of the qubit Bohr frequencies. 



than the qubit-qubit coherent dynamics, since at 20 mK 
it almost disappears. Anyway its presence can be still de- 
tected at temperatures of the order of 10 mK, which are 
reasonable values in the case of two SQUIDs [TJHE]. 

In figure 4 we plot the stationary concurrence C(oo) 
as a function of both T and A. While at zero temperature 



4 Discussion and concluding remarks 

We have presented a general solution of a decay model 
for two interacting qubits, each one coupled to a different 
bosonic reservoir, previously solved in the zero-temperature 
only [35] . The analysis reported here clearly shows that the 
counter-rotating terms in the qubit-qubit interaction, usu- 
ally neglected via a rotating wave approximation, cause 




Fig. 5. Stationary concurrence C(oo) as a function of the bath 
temperatures Ti and T2 for ai = «2 = 10 _3 wi. 
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Fig. 6. Contour plots of the stationary concurrence C(oo) as 
a function of the bath temperatures Ti and Ti when (a) ol\ = 
a.2 = 10~ 3 oji and (b) «i = 10 «2 = 10~ 2 a>i. The white color 



corresponds to C(oo) 
to C(oo) = 0. 



0.45, while the black color corresponds 



two distinct physical effects, still observable at reason- 
able experimental temperatures (typically of the order of 
10 -7- 20 mK). First, qubit-qubit coherent oscillations can 
appear even when the system starts in a superposition of 
the states |00) and |11), a case wherein there would be no 
oscillations in absence of counter-rotating terms. Second, 
and more important, the counter- rotating terms give rise 
to stationary entanglement, which is independent from the 
initial state, a fact that can be exploited to get robust and 
long-lasting entanglement on demand. 

An important point in our analysis is the role of the 
qubit-qubit coupling strength A in the creation of the sta- 
tionary entanglement. Contrarily to what intuition sug- 
gests, it is not useful to increase A indefinitely: after a 
certain value, the stationary entanglement created is very 
fragile against thermal noise. Therefore, one can single 
out a range of values of A within which one can get a rea- 
sonable amount of entanglement even with increasing 



temperatures. The extraction of this robust entanglement 
could be useful for quantum information protocols. This 
point, as well as a non-Markovian extension of the theory 
presented here, will be the subject of our future research. 
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